% 读取DICOM文件夹中的数据
dcmFile = './3#.dcm'; % 指定 DICOM 文件路径
info = dicominfo(dcmFile); % 获取DICOM文件信息
data = dicomread(dcmFile); % 读取DICOM数据

% 确定裁剪区域的坐标
x1 = 766; y1 = 203; % 左上角
x2 = 1090; y2 = 368; % 右下角

% 确定总帧数
numFrames = size(data, 4);

% 初始化存储截取数据的变量
allData1 = [];

% 循环处理每一帧
for i = 1:numFrames
    frame = data(:, :, :, i); % 提取当前帧
    croppedFrame = frame(y1:y2, x1:x2); % 截取定义的区域
    allData1(:, :, i) = croppedFrame; % 保存到数组中
end

% 保存为MAT文件
save('3.mat', 'allData1');
disp('数据已保存为3.mat');
